Luttinger liquid theory of purple bronze Li^gMoeOn in the charge regime. 
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Molybdenum purple bronze Lio.gMo60i7 is an exceptional material known to exhibit one dimen- 
sional (ID) properties for energies down to a few meV. This fact seems to be well established both 
in experiments and in band structure theory. We use the unusual, very 1-dimensional band disper- 
sion obtained in ab-initio DFT-LMTO band calculations as our starting point to study the physics 
emerging below 300meV. A dispersion perpendicular to the main dispersive direction is obtained and 
investigated in detail. Based on this, we derive an effective low energy theory within the Tomonaga 
Luttinger liquid (TLL) framework. We estimate the strength of the possible interactions and from 
this deduce the values of the TLL parameters for charge modes. Finally we investigate possible 
instabilities of TLL by deriving renormalization group (RG) equations which allow us to predict the 
size of potential gaps in the spectrum. While 2kF instabilities strongly suppress each other, the AkF 
instabilities cooperate, which paves the way for a possible CDW at the lowest energies. The aim 
of this work is to understand the experimental findings, in particular the ones which are certainly 
lying within the ID regime. We discuss the validity of our ID approach and further perspectives 
for the lower energy phases. 

PACS numbers: 



I. INTRODUCTION. 



The molybdenum purple bronze, LiQ ^Mo^Oyj , is 
a subject of intensive experimental studies already for 
more than two decades [l|, but its unusual properties 
remain unclear. Several very different experimental 
probes have been used: angle resolved photo-emission 
(ARPES)[2], scanning tunnelling microscopy (STM)|3], 
DC [4 and magneto- resistivity Q, thermal conductiv- 
ity [6 1, optical conductivity Q, Nernst signal [3, muons 
spectroscopy d|, X-rays[iol|. thermal expansionjll| , neu- 
tron scattering 12]. Although the main effort of those 
investigations was focused on the nature of a mysterious 
phase transition at around 25K, interesting knowledge 
about higher energy phase was also gathered. Certain 
properties of the one dimensional (ID) metal, the Lut- 
tinger liquid (TLL), have been invoked to explain the 
measured data@,HG3ll4l an d nowadays the presence of 
the ID physics is well established experimentally [TEl fl(| . 
at least in some energy range. 

On the theoretical side, band structure calculations 
have shown a quasi-lD character of molybdenum pur- 
ple bronze. A remarkably simple band structure emerges 
from rather complex crystal structure. At the Fermi sur- 
face there are only two bands, lying very close to each 
other, in the form of flat sheets dispersing well only along 
the b-axis. This gives a hope that purple bronze can in- 
deed be a rare realization of the ID physics. 

The key problem is that several possible mechanisms 
has been invoked to explain the observed properties, 
which made the subject quite unclear and controversial. 
In our opinion the reason for this situation is that each of 
previous attempts was focused only on one out of many 
peculiar properties of Lio^Mo^On and most of them 
searched for an explanation in the low energy regime (be- 
low 5meV), where indeed the properties of Li . 9 Mo 6 Ou 



are the most spectacular. By now, not enough attention 
has been paid even to the parameters of the ID state. 
It is only agreed that it emerges at energies as high as 
250meV. The values of these parameters are the first is- 
sue one must determine before pursuing research towards 
low energy regimes. The unusual physics observed in 
Lio gMogOir at the energy scales of order lOmeV is ob- 
viously a motivation for revising the question of the ID 
physics. In order to establish a proper low energy ef- 
fective theory one has to begin at highest energies and 
the step by step move towards the physics taking place 
around Fermi energy. The first step is to link the results 
of the DFT calculations with the well defined field theory 
describing the experimental results at around 20meV. It 
is this "high energy" regime which must be well under- 
stood first. This is the main task of this paper. 

The plan of this work is as follows. In SecflT] we be- 
gin with a brief introduction of the band structure. In 
the following section Sec lIII Al we propose a tight bind- 
ing model which is able to approximate bands around 
the Fermi energy Ep, however in addition it contains 
also the strong correlation terms, beyond the LDA-DFT. 
In Sec lIII Bl we give basic notions of ID physics, used 
in the rest of the paper. Section SecJlV] is dedicated to 
intra-chain physics: we give values parameterizing strong 
correlations (Sec lIV A"]) , estimate TLL parameters which 
this implies (Sec lIV Bl) and give energy scales for the spin 
sector (Sec lIV D"]) . In section Sec|V]we introduce the inter- 
chain physics. Once again, first we estimate the strength 
of these interactions (Sec IV At and then (Sec lVB|) cast 
as many of them as possible into effective LL descrip- 
tion, now within the ladder framework. Later in Sec lVIl 
we study how the non-linear interaction terms will affect 
the Luttinger liquid parameters and what instabilities 
they can potentially produce. Finally in Sec I VIII we dis- 
cuss our results for LL parameters in a context of the 
experimental findings (Sec lVIl Al . as well as the validity 
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of the ID approach itself in Sec I VII 51 We also discuss in 
Sec lVIIll the role of substitutional disorder in our model. 
The conclusions in Sec lIXI close the paper. 



II. BAND STRUCTURE 

The lattice space group and atomic positions within 
the structure are known from experiment done by On- 
oda et al [10| recently confirmed by the neutron diffrac- 
tion experiment [12]. This rather complicated structure, 
which consist of well separated slabs parallel to b-c plane, 
is presented on FigQJi). Based on this crystallography 
knowledge the electronic structure of L12M012O34 can 
be calculated using density functional theory (DFT) e.g. 
in the local density approximation (LDA). Recent cal- 
culations [I?} lead to results which are globally consis- 
tent with several previous computations [l|| Gil- Over- 
all the band dispersions agree well with the measured 
results obtained by angle resolved photoemission spec- 
troscopy (ARPES) [20]. For example it shows a flatten- 
ing of the two dispersive bands at about 0.4-0.5 eV below 
Ep. Other bands are found at least 0.25 eV below Ep. 
The only visible discrepancies are the relative vertical 
shifts of the bands of order O.leV. The computed ratio 
between the Fermi velocities along b- and c-crystal axis, 
is about 40, compatible with the reported anisotropic ID- 
like resistivity [j, H, 0, Hl[ . The velocity along a is even 
smaller. At low energies (around Fermi energy Ep) all 
calculations gave qualitatively similar results. There are 
only two bands which cross the Fermi energy and they 
are placed very close to each other. They have a strong 
dispersion in the T — Y direction (b-axis) and barely no 
dispersion in the perpendicular directions. 

The two bands originate from a pair of zig-zag chains 
build out of Mo atoms inside O octahedra. The low en- 
ergies (close to Ep) spectral function is located mostly 
on Ad (precisely t2 g ) orbitals of Mo(l) and Mo(4) atoms, 
using notation from FigQJ,) . These two corners of zig-zag 
are symmetrically inequivalent, thus dimerization is pos- 
sible. However their distinction comes more from out of 
chain, than in chain, environment. The difference solely 
within a chain is hard to notice. It is also hard to dis- 
tinguish the corresponding gap on the edge of Brillouin 
zone upon the analysis of the LDA band structure (in 
the following we this assume that it is not larger than 
- O.leV). 

The standard procedure (see Appendix \K§ is to fit 
the DFT band structure with an effective tight-binding 
model. Along the b-axis (along the chains) the LDA 
calculations were done in the reduced Brillouin zone 
(because of the above described presence of two un- 
equivalent Mo sites). Then, in a tight binding approxi- 
mation, we expect that the band will be back-folded with 
a gap at the boundary of a reduced Brillouin zone which 
corresponds to the above mentioned strength of dimer- 
ization (not larger than ~ Q.leV). In Li2Moi2034 case, 
while it is rather easy to distinguish the lower half of 



dispersion bands (however there are some peculiarities 
at k c « 7r), the upper bands hybridizes strongly with 
other bands (originating mostly from oxygen orbitals). 
Their proper identification is then difficult. The band 
gap seems to be quite small, however the presence the 
hybridization makes the estimate quite difficult. In the 
following we assume that it is < 10%i. 

Along the c— axis the dispersion relation is quite un- 
usual (there is a node at zero momenta along this axis) 
and this peculiar feature appears within all independent 
DFT calculations (l7Hl9j. Here we are also dealing with 
the reduced Brillouin zone, but the presence of the node 
means that the standard back-folding does not apply at 
all. One has to try possible combinations of cosines (see 
Appendix |A"| . keeping in mind two facts: the unequal 
distances for intra- and inter-ladder hopping and the 
large inter-ladder distance which implies that the next- 
nearest neighbor hopping must be very small (only the 
light dashed lines on FigH^) could give non-zero contri- 
bution). 

We were able to deduce (see Appendix [S] for details) 
that the hopping in the perpendicular direction tj_ is 
~ 15meV (certainly smaller than 30meV (« 300A')). 
One should also notice a large frustration, change of sign, 
when hopping to next-nearest chain (interladder hop- 
ping) with the amplitude of this hopping ~ lOmeV. This 
change of sign can be ascribed to a phase shift acquired 
when hopping between pairs of chains, between two Mo\ 
atoms via M02 octahedra (see Figfiji ). 

Along the a-axis we have well separated slabs, with a 
void between them filled by Li atoms and Mo(6), Mo (3) 
tetrahedra. Electrons residing there (if any) stay on en- 
ergies few eVs away from Ep. This explains why LDA 
give a very low dispersion in this direction. 

Significant deviations from LDA were seen in ARPES 
only below 0.2eV. Then the two bands seems to merge. 
This is also the energy scale revealed by optics Q: it 
shows a formation of a first plasmon edge for electrons 
along b-axis (T — Y direction) and a gap value in perpen- 
dicular direction. Thus it sounds reasonable to take it as 
a point where the ID physics forms. This is the starting 
point (on the high energy side) of the present study. 



III. EFFECTIVE LOW ENERGY 
HAMILTONIANS 

A. Tight-binding fermionic hamiltonian 

The LDA-DFT results presented in Sec. HI treat 
electron-electron interactions at the mean field level. It 
is assumed that a single electron moves in an average 
electrostatic potential. However it is known that the 
higher order terms in electron-electron interaction can 
bring very different physics, particularly in the case of 
reduced dimensionality. Thus we need to introduce them 
to our model. 

The LDA results shows that a majority of carriers den- 
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FIG. 1: a) The crystal structure of purple bronze, after 
Ref. [l(J. A cut perpendicular to highly conductive b-axis is 
shown. The atoms which belongs to zig-zag chains are indi- 
cated, b) A simplified structure of the low energy model. We 
show the top view on b-c plane. The t and t± hoppings as 
given in the first line of hamiltonian EqfJ] are shown, sites f 
and 4 lie within the ladder, site 2 is outside it but t± hopping 
most likely goes through this site. 



sity at the Fermi energy is localized on zig-zag chains 
formed by the Mo(l) and Mo(4) octahedra (see Fig.[T]})). 
We can safely assume that the low energy physics is de- 
scribed by the dynamics of these electrons. Then the low 



energy hamiltonian formally writes: 

H = -t\ V ct c - -t\ V c t c . - 

1 / j 1,(7 2+0,(7 * / j 1,<T 2 + 0,(7 

i^everib ,cr i£oddb,o 
i£even c ,cr i^odd c .a 

- t a ^2c\ a c i+s . a + h.c.+ 

i,cr 

(1) 

where the vectors b = [0,6,0], c = [0,0, c], a = [a, 0, 0] 
define the Liq^MoqOiq crystal lattice, summation runs 
over all ladder sites and depends on the parity of a given 
site index along directions b and c. The V n LDA(r m ~ r n) 
are the electron-electron interactions not included in the 
(mean-field) DFT calculation and the sum goes through 
positions of all carriers. The t a ~ c are hopping parame- 
ters along each crystal axis, as estimated in the previous 
section. As discussed there t a « (down to O.lmeV) so 
we can neglect it for the "high energy" range we are in- 
terested in. The electrons are moving exclusively within 
one slab (and mostly along the zig-zag chains). In Eq. [T] 
we kept only nearest neighbor hopping, since they are ex- 
pected to be the dominant ones. The t\ 2 indicate a pos- 
sibility of dimerization, namely different hoppings along 
i-th axis for even and odd bonds. Because Mo(l) and 
Mo(4) are crystallographically different we included the 
possibility that we are dealing with a dimerised chain, 
half-filled in the reduced Brillouin zone. We denoted the 
two corresponding hopping by t \ and t\ . As explained in 
the previous section we take t\ « t\ (with 10% accuracy). 
The t\ 2 describe intra- and inter-ladder hopping respec- 
tively (see Appendix E| . To lighten the notations we will 
denote from now on t b by t and t c by t±, to emphasize 
that they correspond respectively to hopping along and 
perpendicular to the chain direction. 

The resulting simplified tight-binding model is shown 
on Fig.[T|3). We simplify the interactions in a similar way, 
by explicitly considering the intra-chain and inter-chain 
parts of the interactions. The microscopic Hamiltonian 
we take is given by 

H = -t e],<r C J>-*l-L 4,<r C Ji,a-*2_L ^ C V C h,° 

(i,j),<7 (i,jl),<* (i,j2),<7 

m ' n J^m.a^n a' ^7i,a' Cm,er~r 
+ ^ ] Voutiym ~ r n) c m,CT C n,CT' Cn . CT ' C ™ l > CT (2) 

where t = 0.8eV (see Table JTTT|) is hopping between near- 
est neighboring i,j (denoted molybdenum atoms 
along zig-zag chains (b-axis) and ii,2_L ar e a hoppings 
in a perpendicular direction, within a slab (c-axis), be- 
tween nearest-neighbors, I is intra-ladder while 2 is inter- 
ladder. The values of t± t 2± are estimated in Appendix 
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|A1 and summarized in Tab lllll As already mentioned 
t± < 20meV, so we can treat it as a perturbation on 
the top of the intra-chain physics. 

The strong correlations (last two lines in Eqj2]) arc 
usually parameterized by several quantities, which en- 
ter into the effective hamiltonian: U is the local on- 
site interaction between charge densities of opposite spin 
rii-f- = c^Cj-f (Hubbard term), Vi n (r) is the interaction of 
two carriers placed inside the same chain at a distance 
r ^ 0, V out (r) is the interaction of two carriers placed 
in two different chains in a distance r (both Vi n (r) and 
V ou t(r) are defined assuming an environment for which 
the LDA screening is included). We will discuss their 
strength in Sec. IIV Al 



B. Bosonic field theory of a ID system 

In the following sections, it will be shown that the in- 
terchain couplings (both hopping and interactions) are 
reasonably weak compared to the intra-chain ones. We 
can thus anticipate the need to describe a strongly inter- 
acting one dimensional systems. In such a case, a very 
convenient formalism to incorporate the strong intra- 
chain interactions from the start is provided by the Lut- 
tinger liquid formalism (22j. This formalism allows to de- 
scribe non- Fermi liquid state (Tomonaga-Luttinger liquid 
(TLL)) that occurs in one dimension as a result of the 
interactions. 

The low energy dynamics of such a state is captured 
by collective, bosonic modes representing charge and spin 
density fluctuations. These fluctuations are connected to 
the two fields (j> v (x), where v = p for charge fluctuations 
and v = a for spin fluctuations (linked with the fluctua- 
tions of charge and spin density) . These fields have two 
canonically conjugate fields 9 v {x) linked with the respec- 
tive currents fluctuations. In terms of these collective 
variables the hamiltonian reads 22] : 

H ° = E/ ^[(u v K v )^n v f + (^-)(d x cj> v f\ (3) 

where IT„ (x) — -V9 v (x). All the intra-chain interac- 
tions conserving the momentum are now only fixing the 
precise values of the parameters u v (the velocities of the 
corresponding charge or spin density modes) or K v (the 
Luttinger parameters which control the decay of the var- 
ious correlation functions). 

In particular the local spectral function behaves as a 
power law A(x — 0,u,T = 0) ~w° with the following 
(interaction dependent) exponent: 

a = 1 (4) 

For spin-rotationally invariant case K a = 1, thus inter- 
actions, leading to the anomalous behavior, affects only 
the charge sector. 



We can thus take as our starting point for the chain 
physics, and study in the bosonized representation the 
effects of the various inter-chain coupling terms. This 
will ensure that at least the upper energy scales of the 
"high energy" regime will be properly treated. The first 
step to determine the Luttinger parameters u u and K v 
is to estimate the strength of the various interactions, 
which we do in the following section (Sec lIV AJ) . At ener- 
gies lower than the Fermi energy Ep only two scattering 
channels are allowed q ~ and q ss 2k f (plus eventu- 
ally higher harmonics). All density-density interactions, 
with small q momentum exchange, can be incorporated 
into the hamiltonian given by Eqf5] They contribute to a 
highly non-trivial dependence of the TLL parameter K p 
(see Sec UVBl) . 

The remaining interaction terms produce non-solvable 
hamiltonians of Sine-Gordon type, which in general is a 
functional F[) of cosine terms: 

H cos = F[cos{V8pcl) u (x)),cos(VSq9 u (x))} (5) 

where p, q indicate higher harmonics (scattering with 
larger momenta exchange). Due to these terms the total 
problem allows only for approximate solution (at least 
in terms of continuous fields theory) . The terms present 
in H cos are derived from large momentum exchange in- 
teractions (of both intra- and inter- chain origin). This 
is why their presence is analyzed carefully in the further 
sections. They are usually treated using Renormalization 
Group (RG) transformations, which allows to extract the 
terms that affect the most the TLL physics. This usually 
enables one to find the existence of gaps in the spectrum 
of bosonic modes, sometimes even to predict the ground 
state phase diagram. This kind of approach will be ap- 
plied in SeclVTl 

IV. INTRA-CHAIN PHYSICS 

Let us first consider the physics taking place within a 
single chain. We want to obtain (in Sec lIV A]) the values 
of the various intra-chain interactions in ([2]) and from 
there (in Sec lIV B[) the values of the Luttinger parameters 
for the Hamiltonian (J3J. 

A. Strength of interactions 

1. U: local on-site interaction 

A fully self-consistent calculation has been performed 
in [lj| and a value U = 6.4eV was found. This value 
seems at first sight surprisingly large (for Ad electrons), 
but a recent study using constrained-RPA [23| allows to 
understand this result. Although for bulk Mo U s» 3.86^, 
it was convincingly shown that the suppression of the 
pure Coulomb value Uq m lAeV is mostly due to effi- 
cient screening in 3D of the d-electrons In our case (as 
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discussed in details below) the system is underscreened, 
which entails also a significantly reduced plasmon fre- 
quency in comparison with the pure Mo (from ~ 15eV 
to 0.65eV [23]). Thus the value U = 6AeV is justified 
or probably even modest. Then, if we take previously 
estimated hopping t, we get U/t s» 8. 

Note that this makes the local repulsion by far the 
largest energy scale in the problem. 



2. Vin-' interaction inside a chain 

In hamiltonian Eq[2] we defined intra-chain non-local 
interaction Vj n (r). In this subsection we will discuss its 
strength both in real and in momentum space Vi n (q). We 
are going to use shorthand notation V, n = Vi n (r = 6/2), 
as this is usually the single parameter which enters into 
so called U — V model for a ID chain. This parameter is 
much more difficult to estimate than V ou t (see Sec IV A[) 
as it involves the dynamics of the ID metal. Rather than 
trying to estimate it directly from the interaction itself, 
we will simply adjust the value to be put in © and (|3J) 
in order to reproduce the experimental data on optical 
spectroscopy. 

The idea[24|,[25[ is based on two independent estimates 
of kinetic energy in the system. The first one is related 
(by the optical sum rule) to the plasmon (edge) frequency 
I p = ujp and gives us the total possible kinetic energy 
available for the given number of carriers. The same 
quantity can be computed as an integral I a of the optical 
conductivity <r(w) (taken only over the highest conduct- 
ing band). The point is that the second estimate gives 
us the real kinetic energy, renormalized by interactions. 

Fortunately the necessary data is available for 
Lio.gMoO. the value of of the plasmon frequency can 
be read out from Fig. 2 in R ef-B and gives uj p = 0.65eV 
(w p extracted from the LDA [17] is even larger pa l.leV). 
The sum rule integral I a was in fact already evaluated in 
Fig. 3 of Ref. Q • Determining the value of the first sat- 
uration from the plot has of course some error attached 
to it. We take I a ps 0.35 (±5%). This value is rea- 
sonable because at higher energies one expects that the 
other bands start to intervene. The rest of the proce- 
dure is straightforward, the ratio of bare and interaction 
suppressed kinetic energy I a /I p equals to ~ 0.83. Then 
the results of [25| , together with the previously estimated 
value of U allows to predict V, n /t = 1.2 which means 
V in = 0-95eV (± > 5%). 

In addition to the interactions discussed up to now 
(which fit within the U—V model for a single chain) there 
can also exist interactions between more distant neigh- 
bors (see Fig. [2]). These interactions are non-zero be- 
cause of the under-screening, a characteristic property of 
the ID systems. When the separation is larger than that 
between nearest neighbor we can approximate the inter- 
action by a continuum limit Vi n (r > 6) = VcouiMAM 
where 6 = 5.523A is a lattice constant along the b- 
axis[I3| and e(r) is a real space dielectric constant. In 



2D and 3D systems the electrostatic potential is usually 
well screened thus the Thomas-Fermi approximation, al- 
ready used in the DFT calculation, is sufficient. In our 
problem this means that for the very long length scales 
r 3> a the potential is well screened, V(r 3> a) = 
(where a — 12.762A is the largest lattice constant, along 
the a-axis, perpendicular to the slabs planes), because 
then the network of the zig-zag chains can be thought as 
a homogenous 3D bulk. For shorter distances between 
the interacting electrons r ps a and then the distance be- 
tween the neighboring chains does matter. The relation 
a = 2.16 implies that there are a few non-zero interaction 
terms Vi n {r) for a > r > 6. 

The screened potential inside ID wire behaves approx- 
imately like ln(l/r) [22j, same relation holds for the 
screening induced by the presence of other wires [26| . 
In total we approximate that Vi n {q — 0) will get renor- 
malized by a factor (1 + ln(a/6)) ps 1.6. This gives a 
reasonable estimate for Vi n (q — 0), however one has to 
be aware that for the value of charge mode TLL param- 
eter the whole q structure has its importance. 

The values of real and momentum space interactions 
are summarized in tables Tab|T]and TabHTl 



B. Values of TLL parameter in the charge sector 

Clearly the density-density interactions caused by U 
and V are not small perturbations compared to the ki- 
netic energy. In fact they are the largest energy scale 
in the problem, larger than hopping along the chains. 
Thus they shall strongly affect the charge sector, they are 
strong enough to give rise to well pronounced non-Fermi 
liquid properties. We thus need to estimate the TLL pa- 
rameters entering ([3]). As we will prove in later sections 
the following hierarchy of energies holds U 3> Vi n > t > 
V ou t > Jeff (where J e // is an effective superexchange 
which determines energy scales for spin sector). 

We thus start by solving the TLL problem for a single 
chain and will consider later the interchain couplings. 

Even when U ^> Vi„ > t we stay [13] in a gapless 
phase within the TLL universality class, but computing 
the values of the Luttinger liquid parameters K p beyond 
the weak coupling limit is usually a very difficult problem 
[22j |. since one cannot compute them perturbatively. 

In the present case we have to deal with a quarter 
filled system with a very small doping (S around 2%), 
and a small dimerization (as we mentioned in Sec lHI Al 
it cannot be larger than 10% we take moderately 5% in 
any future calculations) and several non-zero interaction 
terms Vi n (r > 6/2). As a first approximation we decided 
to analyze the physics of these zig-zag chains by looking 
at a quarter filled extended Hubbard chain with a local 
interaction U and a nearest neighbor one V. 

If U is the larges energy scale in the problem then 
the charge sector can be mapped onto a spinless chain 
of fermions (or an XXZ spin ch ain) for which the exact 
TLL parameters are known [22l [28] . For a certain value 
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take K p = 1/3 for further calculations. This value corre- 
sponds to the case for which the 4k p charge fluctuations 
decay with the same exponent than the 2k p charge and 
spin density fluctuations [22[. 



where Vf n {U) is the critical value of intra-chain nearest 
neighbor interaction for a given value U. When U — > 
oo the mapping is exact and the critical value is known 
VS,(oo) — 2t. In this case we get K p KXZ = 0.3 which sets 
the lower limit for the TLL parameter. This is a rough 
estimate, but suggests that for further analysis we should 
use approaches which are valid for the gapless phase. 

A more precise estimate can be obtained from a nu- 
merical study of the extended Hubbard model, and an 
estimation of the TLL parameters in the usual way via 
thermodynamic quantities (28l - l30j . From the relevant 
plots one reads out that our model is located somewhere 
in the range K p £ (0.3,0.37). One can try to give an 
even better estimate using EqJS] even when U is not in- 
finite (but still with U much larger that any other en- 
ergy scale in the problem). It is commonly believed that 
the value V£ n (U) extracted from the numerics can give 
quite good approximation when substituted into ([5]) • Es- 
timates of critical Vf n (U) are usually done with higher 
precision, than for arbitrary U, V . In our case U = 8t 
implies V&(U) = 2M [H, Eland V&(U) = 2.75t from 
older work by Mila and Zotos[29]. This gives K p = 0.330 
(with V/V c = 0.72) or K p = 0.340. We have also fol- 
lowed the critical scaling analysis proposed in [3lj and 
found similar value K p = 0.328. 

The exact solutions e.g. Bethe ansatz are available for 
a few special cases. If we include a (quite weak) dimer- 
ization and (very small) doping then the problem is lo- 
cated far away from any integrable model. However the 
influence of both these perturbations are known. The 
dimerization is able to lower a bit value of K p , this ef- 
fect can be present in particular in our case when we are 
not far from V? n [3(| ■ The doping has the opposite ef- 
fect, however when one is not very close to critical point 
K c = 1/4 and close to commensurate case then the effect 
is negligible. 

To conclude: the knowledge about the U, Vi n we have 
collected above allowed us to give a several estimates 
based on complementary numerical calculations of U — V 
models. All these estimates gives K p £ (0.3; 0.36) (we al- 
lowed for 10% error in the U, Vin due to uncertainty of 
parameters). This is also in agreement with a very re- 
cent study for a system with finite range interaction [32j: 
K p « 1/3 when we take similar value of on-site term 
and account for the presence of Vi n (r) interactions up to 
four-five Mo sites. However in this last paper interactions 
with an exponential character were studied, which means 
that they decay much faster than those in our problem. 

Note that this K p value is reasonably close to the one 
K* = 1/4 that would lead to a quarter- filled Mott in- 
sulator in the presence of an infinitesimal Vi n . The sys- 
tem will be thus very susceptible to the precise value 
of Vi n . Given the accuracy of our estimation we will 



C. The remaining interaction terms 

By now we have studied interactions in real space 
(which we believe can be useful for numerical studies), 
while cosine terms beyond the TLL general expression Q 
are defined in the reciprocal space representation. These 
values of scattering for large momenta exchange will be 
used as an input for the discussion done in Sec. I VII 

The Hubbard U interaction has a form of a delta func- 
tion in real space thus in momentum space it contributes 
equally to small and large momenta exchange scattering 
(provided they are the intra-chain ones). The previous 
considerations also imply that, the Fourier transformed 
Vi n {q) is a weakly decaying function (slower than loga- 
rithm) , so it affects the amplitude of scattering processes 
with large momenta exchange. As it was already ex- 
pressed above the intrachain interactions are much larger 
then t thus it is not straightforward to obtain the value 
of backscattering they cause. In fact the K p — 1/3 pa- 
rameter given above, known from numerical studies of 
similar models, is a fixed point value K* which means 
that it contains a part of these contributions. Precisely 
it is the part included in the simplified U — V model. The 
longer range, slower decaying interaction leads to smaller- 
value of K* . This suggests that the proper value of K* 
in our problem which contains long range interaction can 
be even smaller. 

In our case the charge sector is renormalized by the 
possible umklapp terms 53, so the last statement can be 
rephrased: the smaller final value of K* can be linked 
with larger value of initial, bare g®. The precise esti- 
mation of the amplitude of g® is a very difficult task, the 
detail discussion of the bare large-q, intra-chain, term will 
be done in Sec. I VII 

In addition, the e(r) when computed using RPA (be- 
yond standard DFT) gives a well known phenomena, the 
Friedel oscillations (in real space). The effect comes from 
the peculiar screening (singular susceptibility) at large 
momenta q = 2kp. In our case, due to the value of kp 
along chains, it gives an extra gain of energy for electrons 
located every second site. This affects the large momen- 
tum exchange part of interaction Vi n (q ~ Akp). The 
value of this gain can be estimated using the fact that, 
for 3D metal, the additional oscillating part V(r) ~ r~ 3 
thus we give an estimation Vpna < V tn {r = b)/8 = O.leV. 
It is not as large as U or V, but can be significant if we 
compare it with t±. 

The U-V model that we have considered so far (in par- 
ticular in Scc lIVBj) is of course only an approximation 
of the intra-chain physics and interactions. It should 
capture most of the effects at the energy scales we are 
considering, however it may slightly underestimate the 
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strength of interactions. 



D. Spin mode 

All the interactions considered up to now were con- 
nected with the charge sector. In ID systems, because of 
the spin-charge separation, electrons spin degree of free- 
dom should be discussed separately. For a half filled chain 
the knowledge about t and U allows to estimate the spin- 
spin exchange (superexchange) constant J = tjj. This 
determines the energies at which spin sector starts to play 
a role. The problem of purple bronze is more complex 
since the compound is quarter filled. 

For such cases the formula for superexchange interac- 
tion can be still obtained from second order perturbation 
theory goj. It reads: 



J eff ~ 



4*| 



8ti + 2U + V in - 2y/(U ~V m ) 2 + W{ 



(7) 



where we have taken into a account the fact that due to 
dimerization there are two slightly different alternating 
hopping t\ and *2 along the chain taking *2 — *i ~ 0.05* 
and using the values of U and Vi n obtained in Sec II V Al 
we obtain: 



Jeff = 0.2eV 



(8) 



It is indeed significantly smaller than * (and also the 
charge sector interactions U,Vi n ) which implies that 
charge dynamics will dominate the 0.2 — 0.02el^ energy 
range. However it is a non-negligible value in the sense 
that even for energies as high as ui ~ O.leV (in the middle 
of the considered "high energy" regime) the spin excita- 
tions are coherent and dispersion is linear thus a TLL de- 
scription with spin and charge modes is applicable. The 
spin- incoherent TLL 33] is not an appropriate framework 
for our problem. 

For a spin-rotational, SU(2) invariant model we have 
K a (T = 0) = 1. For the interacting case the spin ve- 
locity u a is smaller than vp ~ 2*, in particular v a (U — > 
oo) — > J. This last value is comparable with the one 
observed in experiment [131] . The only candidate to break 
the SU(2) symmetry would be spin-orbit coupling Dls 
on the heaviest atom, molybdenum. A series of experi- 
mental and LDA studies allow to set its value in bulk bec 
Mo[13] to Dls = lOOmeV. This value was obtained for 
the A point of a Brillouin zone in the Mo-bcc crystal. In 
our problem it should be smaller because the active elec- 
trons have mostly ti g character (with a larger /-number 
value). Thus Dls can be treated as a perturbation for 
J, contrary to U for t in a charge sector. This implies 
that K a deviates from the non-interacting value K a = 1 
much less than K p . To be precise from weak coupling 
theory (applicable in this case) we know(22j: 



Dls 
2nJ 



1 - 



Dls 
2ttJ 



(9) 



from which we predict K a =1.1 which can increase the 
Green's function exponent a only by 0.01. The spin sec- 
tor thus cannot be responsible for the experimentally ob- 
served values of the a exponent (defined in Eq llip which 
are of order ss 0.5. However the influence of the spin-orbit 
coupling should be taken into account when one studies 
the physics taking place lower energy scales, which are 
beyond the scope of this paper. 



V. PHYSICS BETWEEN THE CHAINS 

After giving the description of a single chain physics we 
move beyond this model and study the strength (and the 
role) of interactions between carriers moving in different 
chains. 



Vout'- the values of interactions between the 
chains 



The density-density interactions can be estimated as a 
standard Coulomb interaction between charges in a 3D 
dielectric, as it was done in Ref.[l8j]. Those authors gave 
estimate for a value of interaction from standard electro- 
static Coulomb law (with the simplest static screening): 



V out (r = a) 



1 



KEQa 



(10) 



which in fact gives the value for two electrons in two 
different 2D slabs. We are more interested in interac- 
tions inside the slab which is obviously larger because 
(the smallest) interchain distance is then c/2 = 0.4a. A 
more fundamental problem is the value of the dielectric 
permittivity k. In the previous work [18| the bulk R sa 10 
was used, which is typical for bulk semiconductors with 
similar value of a gap for oxygen states(Ao ~ 2eV). Tak- 
ing into account the very weak metallic character along 
the c-axis the semiconductor approximation is correct. 
For further neighbors, with many oxygen atoms in be- 
tween, one can take the bulk R value. However for the 
nearest chains, as there is only one raw of oxygen atoms in 
the space between them, such large n value overestimates 
the screening provided by sparse environment [43]. It is 
then convenient to take a function such that k{c/2) 
is reduced by 30% in comparison with R. For larger r dis- 
tances at approximately r « a it saturates to bulk value 
R. With this set of values we estimate V ou t — Q.hheV . 
For r — > oo the metallic character along b-axis intervene 
and k(t) —> oo. Thus the large distance interaction is 
strongly suppressed. In such a case the above estimated 
value for Vout can be taken as the density-density interac- 
tion (V ou t(q ~ 0)) between the nearest chains. The above 
estimate was done between two nearest chains which form 
a pair (see Fig. [2]). On the other side of each chain there is 
another neighbor which is placed two times further. Ac- 
cording to (fT0|) these interactions V ou t2 (we defined the 
inter-ladder term in analogy with e.g. t C 2 in Eq |Al[) are 
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FIG. 2: (color online) A real space image of strong inter- 
actions terms possible in purple bronze; the top view of a 
conducting slab is given. Interactions between an electron lo- 
cated on a dot-site and charges on different sites are indicated 
as red dashed lines of different thicknesses. Both Vi„ (dashed 
lines within a black chain) and V ou t (dashed lines between the 
black and blue chains) are shown. The numbering of sites is 
like on FigfT] 



at least two times smaller. We then estimate (keeping 
n(r) in mind) V out 2 — 0.2eV. 

The treatment of terms with large momentum ex- 
change is more complex. There is a q — 2kp term 
which locks the interchain density wave, the so called 
ir — CDW, which is shown on the Fig. ■ From the fig- 
ure it is clear that the distances r„ = \J (c/2) 2 + (2/cf) -2 
between charges in such configuration are rather large 
(around three times larger than the distance entering to 
previous calculation for V ou t(q = 0)) thus the screening 
is quite efficient and n[r v ] is definitely not reduced, but 
probably even enhanced with respect to the bulk value. 
In fact not only further chains, but also further slabs can 
intervene (because r n > a). 

There will be either a very efficient screening (like in 
a metal), or Coulomb potential approximation (but with 
enhanced k) is still applicable. We assume, optimisti- 
cally, the second case and use again (fl"0|) . Two ways of 
proceeding are possible when one is interested in the stag- 
gered component of the interaction between two chains. 
In real space (direct application of (jTUJ)) one computes in- 
teractions with a linear set of dipoles. Due to increase of 
n(r) the interaction (with further dipoles) decays rapidly, 
which makes this straightforward approach quite tedious. 
In the reciprocal space approach one must take into ac- 
count the fact that Fourier transform of V ou t (q) does de- 
cay in momentum space (it will be a 1/q decay in a 2D 
case corresponding to a separate slab). The second ap- 
proach is simpler when one notice that the previously 
computed V out — 0.55eV correspond to q = ± (a/2) -1 
(a distance between slabs sets the large distance cut-off). 
Then the 1/q scaling allows to estimate that q = 2kp 



a) 




FIG. 3: A schematic side view of a ladder structure (which 
enters to the low energy description of purple bronze) with: 
a) the 2kF it— DW instability shown (The light, dashed line 
correspond to the distance r n ) b) the 4&f n— DW instability 
shown. 

term will be suppressed by an extra factor four. Overall 
we have: 

• a factor 1.5 -j- 2 from the value ft(rv) 

• a factor 1.5 -j- 3 from the distance r„ 

• a factor s» 4 because we compute the staggered 
(large q) component 

Taking all these factors into account we estimate that 2k p 
component is reduced by at least one order of magnitude 
which means V out {q = 2kp) < 0.05eV. 

There is also the q = 4k p term. Now the distance 
is smaller, thus k is unaffected. However this term is 
strongly suppressed due to above mentioned Coulomb 
character of V out (q)- At short distances corresponding to 
r = (4k f) _1 the interaction may even have a 3D charac- 
ter (with no screening) thus V out ~ q 2 . By a similar rea- 
soning like above we get V ou t(q = 4kp) m 0.05eV. The 
smallness of inter-chain large momenta terms (in com- 
parison with density-density term) implies that the local 
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TABLE I: The real space values of strong coupling parameters 



u 


V in (r = b/2) 


V in (r = b) 


V in {r = 2b) 


V out (r = c) 


spr 

V out2 {r = g^. 


6.4 


0.95(± > 5%) 


0.5 


0.2 


0.55 


— — aefr 



The respective values are: K p s ~ 0.30 and iC^ w 0.38. 
For energies well below V ou t the hamiltonian of the charge 
;or is given as a sum of two bosonic modes plus a 
-linear part H coa caused by large-q (exchange) inter- 
ions: 



TABLE II: The reciprocal space values of strong coupling 
parameters (given in eV). 



H, 



eff 



E 



Vin{q = o) 


V in {q = 2k F ) 


V out (q = 0) 


V out {q = 2k F ) 


V out (q = 4k F ) 


1.6 


see SeclVI Bl 


0.55 


< 0.05 


0.05 



field corrections are small, so our mean field model of 
screening is self-consistent. The influence of the inter- 
chain exchange interactions (both 2kp and 4fcp) on ID 
physics will be discussed in Sec. IVII 

To summarize the results for values of the strong cou- 
pling parameters found in the last two section we present 
all parameters in the tables Table U and Tabic [TTJ 



B. Luttinger liquid framework 

The density-density interaction between the two neigh- 
boring chains V ou t(q = 0) can be included on non- 
perturbative level. Because of the hierarchy of energies 
the presence of V out should be included on the top of ID 
chain K p . The zig-zag chains are grouped in pairs as 
shown on Fig. [T|d) . Inside each pair we have two short 
links with interaction V ou t while in between them there 
is only one link with much smaller interaction V ou t2 ■ The 
ladder picture is then justified. 

In the case of the ladder-like system (which means 
that there is either much stronger intra-ladder hopping 
tix 3> f2_L or interaction) two more modes must be in- 
troduced. The problem can then be expressed in two 
different possible basis. One may work either in the 
chain basis with v — oT,cr2 spin modes or in the to- 
tal/transverse basis with symmetric and antisymmetric 
spin modes v = aS, a A (and respectively the same for 
the charge sector). In general the Greens function a ex- 
ponent can be expressed as: 



2N 



(11) 



where N is a number of modes. A significant V ou t (q = 
0) gives a preference for symmetric and antisymmetric 
modes as it can be diagonalized in this basis. 

Thus within the ladder picture the two degenerate 
charge modes (in the two chains of the ladder) split into 
two holons with two different velocities. One of them 
(symmetric) describes the fluctuations of the total charge 
density, while the other one (anti-symmetric) describe 
the relative fluctuations in two different zig-zag chains. 
Their values are: 



K. 



K n 



pS,A 



1 =F K p V out /(TTV F ) 



(12) 



pt(uM(?rTl„f + (^)(d x ^) 2 ] 

~T~ Hcc?, 



(13) 



where as explained in Sec lIIIBl now the summation goes 
over four modes v — pS, pA, aS, a A. If we neglect H cos 
then the single-hole propagator (relevant for ARPES), 
the density Greens function with chirality #: 



G<(x,t) = (iPl <7 (x,t)^v(0,0)) 



(14) 



is known [3 5| and can be written as a simple general- 
ization of the standard TLL result, as a product of four 
gaplcss modes. For example for the right going hole: 



where: 



0C, = [{x - v v t)] Q " [{x + v v t)f 



(15) 



with the exponents C„ = (K v + K~ 1 )/8. For a nu- 
merical values of these exponents in our model and their 
implications see Sec lVII M 

In the case when H cos is able to open a gap in a given 
mode, the respective term in the product in Eg IT5l should 
be substituted by a modified Bcssel function of a second 
kind, which gives the expected exponential decay of cor- 
relation function. 

The inter-ladder interaction V ou ti are significantly 
smaller and thus will not bring any novel physics at 
higher energy scales (above lOmeV), so this leave them 
out of the scope of this work. 



VI. RENORMALIZATION GROUP STUDY 

A. Statement of the problem with the inter-chain 
operators 

The interactions with small momentum exchange can 
be absorbed in a definition of TLL parameters, but ob- 
viously there are also scattering channels with the large 
momentum exchange. These generate cosine type terms, 
which in principle can bring us away from TLL univer- 
sality class as defined in EqJ3J There are several terms 
which are expressed as cosine operators in the bosoniza- 
tion language. In addition to previously incorporated 
umklapp terms there are the ones which emerge from 
the inter-chain interactions. It is because usually they 
are the most pertinent for the ladder system. They are: 
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qi = 2kp these are interactions originating from the 
presence of V ou t(q = 2k f)] they have a form 
cos(2cj) pA ) cos(\/20 cr i) cos(v / 20a2) 

qi = 4:kp these are interactions originating from the 
presence of V out {q — 4fcp); they have a form 
g\ cos(40 Pj 4), and (only in the commensurate case) 
the inter-chain umklapp 3" cos(4</> p ,s) 

t± the single particle hopping induce several co- 
sine operators (36j . each of them in the form 

cos(20p j4 )F[cos(v / 2</) (T i), cos(v / 20 (T 2)], where func- 
tional F[] is a linear combination (there are also 
higher order hopping terms, but as they are pro- 
portional to i j_ we can safely neglect them) 

The standard way to treat these terms is deriving, per- 
turbatively (usually on the single- loop level), the renor- 
malization group (RG) equations [22[ . The RG equations 
allows us to predict weather a given term will increase 
and affect the low energy physics (become relevant) as 
the running energy variable A max[uj,T] decreases. 
To be precise the quantity I entering to RG equations 
is defined as A ~ exp(— I) or to be precise I = ln(A/W), 
where W = 2t. This last formula allows to link two quan- 
tities (which we will do throughout this section), however 
one should remember that while I is just a number, A has 
an energy unit. In the following we apply the following 
convention <?J[Ai] = V ou t(q — 4fcp)/A where A = irvp. 

The naive way would be to add RG equations describ- 
ing these terms into the previously determined TLL fixed 
point. Unfortunately, in our particular problem this ap- 
proach is not justified. Our problem can be stated as fol- 
lows: one to incorporate the above given operators into 
the intra-chain, umklapp RG flow, which was already ac- 
counted for. 

We expect that the following physics takes place. 
Around 200meV a gappless TLL appears and later un- 
dergoes the renormalization flow. At this finite en- 
ergy we know the values of interchain terms caused by 
V±_(q = 4fcp) as they were estimated in the previous sec- 
tion, but 33 cannot be taken to be equal to zero (it is not 
yet a fixed point as explained in Scc lVI B|) . If we want to 
treat all instabilities with 4kp periodicity on equal foot- 
ing then we have to begin the flow with non-zero intra- 
chain 33 and see how it will compete or conspire with 
inter-chain instabilities. Estimating the value of the ini- 
tial 33 which we need to substitute into RG equations is a 
highly non-trivial task, our proposition on how to tackle 
this problem is given in appendix [Cj 

B. The intra-chain, umklapp RG flow 

In Sec lIV Bl we have given the values of intra-chain TLL 
parameter K p within U-V model approximation. The 
point is that these are already renormalized, the fixed 
point values K* , and as such already contain influence of 
intra-chain umklapp scattering. 



The fact that this operator is there and affects K* 
becomes clear from the following reasoning. As men- 
tioned above, we are working with weakly dimerized 
chains close to half-filling. In the reduced Brillouin zone 
k F = 0.4877r/& according to DFT and k F = 0.5l7r/& ac- 
cording to ARPES. This discrepancy can be understood. 
One should remember that the precise value of kF can 
vary as it depends on the relative value of chemical po- 
tential inside a chain (determined also by strong corre- 
lations) with respect to the local potential on the Li ion 
(donor of electron) . In any case one must admit that the 
doping is low enough, so in the "high energy regime" the 
zig-zag chain is not able to recognize whether it is doped 
or not. 

In addition an important remark has to be kept in 
mind: the above form of umklapp operator is appropriate 
for a quarter filled chain. This approach holds for weakly 
dimerized zig-zag chain, to be precise for the case when 
amplitude of umklapp interaction is larger than dimer- 
ization. From the analysis of the DFT band structure we 
know that tx — tz < O.leV. After evaluating the strength 
33 we also check if this assumption was consistent. 

In this case the umklapp terms (so called (73), in the 
i-th chain, in the form cos(^/8(/) P i), has to be taken into 
consideration. The RG equation for this instability reads: 

^ = 8 5 3(l/4 -K p ) (16) 

The critical value of TLL parameter for this opera- 
tor is K° — 1/4, so it seems to be irrelevant in our 
problem. However in the case of Berezinskii-Kosterlitz- 
Thouless flow, like the one described by eqUU there ex- 
ists a straight line on the 33 — K p plane, a separatrix 
between relevance and irrelevance regime (gapped and 
critical phase). Essentially this is because the compress- 
ibility of charge mode will also change: 

^=-8glKp (6'ff(l)) (17) 

and for large enough values of initial 33, this equation 
cannot be neglected, the RG flow cannot be assumed to 
be vertical as we did before. Thus the amplitude of bare 
coupling 33 can play a role. As we will see below, in our 
problem the amplitude 33 does matter. 

In the equation Eq |17l we kept the doping dependence 
which is encoded inside the Bessel function of the first 
kind Jq(5(1)). However below we assume that doping is 
essentially equal to zero and we work in a commensurate 
case. This assumption comes from the fact as we work 
in a constant chemical potential (the chain is embedded 
in a whole crystal and one knows the chemical poten- 
tial of such a system from DFT solution) the effective 
doping will get renormalized ~ — g^Jo{5 e ^ {I)) already 
at higher energies. Because of the same reason we ne- 
glect the renormalization of the charge modes' velocity 
~ g^K p J2(8 e * * (I)) ■ From now on all we neglect all dop- 
ing dependence and set 8 e ^ — in the "high energy 
regime" . 



11 



C. The RG analysis of the 4&f terms 

If the inter-chain 4k p terms are considered alone then 
we immediately see that they are less relevant than the 
2kp instabilities described in Sec lVI Dl (it is because 
K pv — 1/2). However there are several reasons why we 
think that the 4kp will be more important and decided 
to investigate it first. The 4k p depend only on charge 
modes, which can be particularly important for energies 
larger than or comparable with J. What is more, as we 
deduced before, K p < 1/3 and then it is the 4/cf CDW 
instabilities which are decaying slower. Finally there is 
an intra-chain umklapp term with the same periodicity 
and large, bare amplitude. 

Section Sec lVI Bl gave us necessary understanding of 
the intra-chain RG flow, now we can proceed and intro- 
duce inter-chain terms. The q 2 = 4kp it— wave scattering 
has the following RG equation: 



dgl 
dl 



2(1 - 2K pA )gl - g 3 gl 



and the inter-chain umklapp: 



dn u 

= 2(1 -2K pS )gl- g 3 gl 



(18) 



(19) 



In the above equations we have already introduced the 
so called mixed terms caused by the presence of g 3 term 
(see below for an explanation). We need to take the 
inter-chain terms together with the standard umklapp 
g 3 to get a complete RG flow. These terms will also 
affect the flow TLL parameters, the K p s,a, the symmet- 
ric/antisymmetric TLL parameters, which are defined in 
different basis that the intra-chain -Kpi, 2 ones (used in 
Eq|16l and II T[) . The intra-chain umklapp scattering may 
be rewritten in this basis using the fact that the com- 
bination of two umklapps can be expressed in a rather 
simple form: 



(20) 



ff 3 [cos(2v / 80 p i) + cos(2a/80 (32 )] = 

2g 3 cos(4(f> pA ) ■ cos(4(f> p s) 

where we assumed that umklapp g 3 is identical in both 
chains 1 and 2. This can be justified by the crystal sym- 
metry argument. 

In the following we work with the perturbation de- 
scribed by Eqf2D]in the symmetric/asymmetric basis. We 
need to re- write the RG equations for the combined intra- 
chain umklapp. Instead of Eg 11 61 and Eg 1 171 now we have: 



dg 3 
dl 



g 3 (2-4(K pS + K pA )) + 2glg] 



dK 



01 



dK, 



pA 



01 



~K 2 pA (8gl + {glf) 



(21) 
(22) 
(23) 



where (in the last two equations) we have distinguished 
the two RG flows of symmetric and antisymmetric K p s 
(instead of single Eq|T7|) and already included their de- 
pendence on inter-chain interactions V ou t- These equa- 
tions should be taken together with Eq[TS] and Eq03)] to 
obtain the full RG flow. Below we show in Figf4]the re- 
sult of a direct integration of this system of differential 
equations as well as a semi-quantitative analysis of RG 
flow between energy scales corresponding to Ai = 0.2eV 
(where TLL is likely to form) and A 2 = 0.02eV. These 
are the limits of interest for the "high energy regime" 
study. The physical reason of this limit will be given 
below. 

One can dispute whether a direct integration is valid 
when g 3 and yi\ (see AppJC]for y» definition) are of order 
0(1). Such large terms could in principle generate sig- 
nificant higher order operators along RG. However there 
are no higher order terms proportional to 173 (n > 2). 
Then the higher orders mixed terms will be either irrel- 
evant (because of large p value in functional Eq(5]) or of 
quite small amplitude. Qualitatively the full RG flow 
can be thought as superposition of two BKT flows, two 
hyperbolas on gi — K v plane. 

Let us first check what are the amplitudes different 
instabilities gi at A2, in the first order approximation. 
If we neglect the flow of K V {1) then the solution of the 
equations Eq[THJ [IT)1|2"T1 is an exponential function: 

ft,[2i*exp(-0] =s„[Ai]exp(»£(-i-ln(Ai/2i))) (24) 

where r\ v is the dimension of the operator <?„, e.g. rf^_ = 
2(1 — 2K* A ) The crucial approximation we made in Eqf2H 
is that we assumed exactly vertical flows and took the 
fixed point values of K*, which tends to overestimate 
the strength of all instabilities. As a result of EqJM] we 
get: 53 [A 2 ]A a O.W, <?I[A 2 ]A a 0.075eV, #[A 3 ]A « 
OMeV. Note that because initially K pS [hi] < K p a[M], 
the inter-chain umklapp <?" renormalizes more than the 
7T— DW one g±, however one has to remember that its 
bare amplitude is smaller due to doping effects. Already 
this simplified reasoning shows that the amplitude of the 
irrelevant term even at A 2 is still larger than amplitudes 
of the relevant terms. This is a peculiarity of our prob- 
lem. The presence of a mixed terms will in fact enhance 
this property: because g\ ~ g\ < g 3 then one can inter- 
pret EqfTBl and EqlTH] as if K c pS A was shifted from 1/2 
downwards. On the other hand in Eqf2T] we see that 
inter-chain instabilities push (weakly) the RG flow of g 3 
towards its separatrix. 

Now let us move to the analysis of the K p s, A RG flows. 
From Eql^U and Eqf53] we immediately realize that intra 
and inter-chain terms support each other in lowering the 
values K p s. A . Taking into account the initial (at Ai) hi- 
erarchy of energies we can be sure that the initial flow 
of K pv ~ g\ will be dominated by g 3 while the reason- 
ing of the previous paragraph showed that around A 2 
both terms are equally important. Initially the flow slows 
down because an irrelevant 173 decreases, but later it can 
speed up again due to relevant inter-chain instabilities 
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FIG. 4: (color online) Three RG flows of K p s (red, three 
bottom curves at I = 0) and K p a (yellow, three top curves 
at I = 0) TLL parameters for different choices of initial pa- 
rameters. The energy scale A2 corresponds to I « 3.5. If all 
K v tend to be constant for increasing Z, then the RG flow is 
so called "vertical". The initial values for three curves are 
(from top to bottom): (0.05, 0.04 , 0.5 , 0.5, 0.4);(0.05, 0.04 
, 0.4 , 0.45, 0.35);(0.06, 0.05 , 0.38 , 0.45, 0.35); where the 
following notations are used (<?1A, g±A, 53A, K p a, K p s). For 
the RG flow which give the bottom K p s,a curves, the insta- 
bilities (not shown) are large gi « 1 at the scale A2 which 
means that the RG flow should be stopped. The other two 
flows do not suffer from this limitation. 



as seen on figure Fig|4] where the result of numerically 
solving of Eqi22knd Eql23l is shown. Estimating quanti- 
tatively -Kps : a [A2] can be achieved in two ways. 

First way (i) assumes the independence of intra- and 
interchain BKT RG flows. We already know the influ- 
ence of 53 term alone, now we want to compute how much 
the K V [A 2 ] would be lowered during the BKT RG flow 
caused only by the inter-chain term, e.g. if we keep only 
g\ term in Eq[23]for K p a\A^. We can use the procedure 
very similar to the one applied in AppiCj just that now we 
are moving towards I — » 00 along the RG trajectory. By 
analogy with AppJC] reasoning we define the flow invari- 
ant A±. In the case of this RG flow g± ss 0.01 <C 1 thus 
A± [y where (y )jj = 0.11 is the distance of the 

bare K° pA = 0.39 to K c pA = 1/2. We then obtain (from an 
analog of Eq |C2[) that t/j| [A2] = 0.15, which means that 
RG flow caused Ayf[A 2 ] = yf[A 2 } - A± m 0.05 of the 
-K/>a[A2] decrease j44[. This additional change of -ftfpA^] 
should be added on the top of the value previously found 
from the U-V model (see Scc llTB]) K* A 1/3. However, 
because of the 53 influence, the A± is not a constant but 
it increases during the flow so the above value of the 
Ayjf[A2] is underestimated. Such a deviation from a sin- 



gle BKT flow, particularly pertinent at lowest energies, 
is clear from hgure Figj5l 

The second way (ii) comes from the fact that, as 
pointed out above, the 33 term dominates most of the 
RG flow between Ai and A 2 . Let us assume that the 
inter-chain terms effectively adds up to the initial ampli- 
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FIG. 5: (color online) A comparison between BKT flow 
caused by a non-zero $3 only (blue, bottom curve) and the 
result of our full RG flow (red, top curve). This is a para- 
metric plot with <j3 [/] on the vertical axis and an average 
(K p s + K p a)/2 on the horizontal axis. The initial values are: 
for blue line (0.0, 0.0 , 0.35 , 0.45, 0.35) and for red (0.06, 
0.05 , 0.35 , 0.45, 0.35) with the same notation like in Figg] 



tude of 53 [Ai] (see EqEUJ): 

gl ff [A x ] « ,g 3 [Ai] + A aux * 2 * V out {4k F )/A (25) 

where the origin of A aux =0.6 factor is explained in the 
end of this section. Taking into account the results of 
AppOwe see that with this new effective amplitude of 
the umklapp the initial point of the RG is located quite 
close to the separatrix of intra-chain flow. The resulting 
value of K* pv is then going to be close to K pv = 1/4. This 
obviously overestimates the change of K pv \A 2 \ because 
effectively in Eq|22]we took a square of a sum instead of 
a sum of squares. There is also a hidden assumption here 
that K p a ~ K p s, but looking at figure Fig|4]we see that 
this works well for all cases. 

From the reasonings (i) (upper limit) and (ii) (lower 
limit) we conclude that within our approximation 
Kps.A S (0.25,0.29). This is in agreement with the re- 
sults presented on Fig |4] and FigJSJ 

Finally we can try to estimate the value of a gap in 
the holon spectrum. We choose to work with g\ as it 
has the strongest tendency to open a gap. We estimate 
the gap in the case if this instability was acting on its 
own when A sa W exp(— I*) where g\ [W exp(— /*)] = 1. 
For the specific cases the I* is known [2 2] for example 
deep inside the gapped phase (self-consistent harmonic 
approximation) or close to the separatrix of RG. Unfor- 
tunately our problem does not belong to any of these, 
because A± ps 0.12 which for A2 energy scale is neither 
very small nor very large. We approximate the flow of 
the considered g\ by BKT flow which can be integrated 
out (see also ApplCj): 



» / 2/1] 

arctan — — — arctan — - 
I A± ) \Aj 



4^- )=Ad 



(26) 
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If we use the fact that yjj « A± and y\\ [2t exp(— I*)] >• A± 
(we are not far from the separatrix) then we get A±l* ~ 
7r/4. In fact this value of I* is quite close to the ones we 
get from numerical integration of the full flow. Then the 
gap is: 



Al = 2tcxp [-7r/(4Aj 



(27) 



which gives ~ 2meV. The low initial value of the 
Vout(q = 4/cf) translates into quite small expected value 
of the gap in the spectrum. 

The lowest energy (below A2) flow of umklapp pro- 
cesses are seriously affected by doping. If the g\ is able 
to open a gap before that then these other RG equations 
will be changed by the presence of the gapped mode. For 
the intra-chain umklapp instability the initial point lies 
very close (esp. if we include terms beyond U-V model) 
to separatrix, which means that it is very susceptible to 
these modifications. 

Usually in such a case one thinks about two competing 
instabilities taking place in two different basis (intra- and 
inter-chain). Our problem is different, the cooperation of 
instabilities is realized. This becomes clear if one assumes 
that he tends to be locked at minimum e.g. correspond- 
ing to V8(/)pi,2 = tt (then cos(v / 8</> p i i 2) = — !)■ For this 
specific value of <jP pl 2 the combined inter-chain terms (see 
Eq |20|) give non-zero and negative value (A aux s» —0.6). 
This corresponds to an additional energy gain caused by 
an auxiliary instability. It is a rare case when the two 
gaps does not exclude but can help each other. In par- 
ticular this validates the assumption made in Eql25l 



D. The analysis of the 2k f terms 

The RG equation for the 2kp backscattering operator 
reads [22}: 



• at lower energies (~ 0.015eU) several terms gen- 
erated by the perpendicular hopping (~ t±) will 
start to intervene. Although they are irrelevant 
(connected with the field 9 p a [36j), still their ini- 
tial amplitude is significantly (three times) larger 
than V ou t(q = 2k p). We then expect that 
physics will be dominated by a competition be- 
tween these two types of cosine terms, as it was 
studied in Ref . [37|, [38[ where the term confinement- 
deconfinement transition was coined to describe the 
physics. 

Based on results from Ref. [13, [Hj], one may expect that 
due to last mechanism the 2k p instability is weakened. 
In any case because K p a <§C K pA — 1 the 2k f— RG flow 
is vertical. This implies that K p a value (and thus also 
Green's function a exponent) is very weakly affected by 
the presence of the 2kp terms. 



VII. DISCUSSION 

A. Comparison with experiment 

In this last part we wish to compare the estimated 
Luttinger liquid parameters with experimental findings. 
Experiments look at low energy (I > 1), long distance be- 
havior (r > fc^ 1 ), which we should keep in mind for the 
rest of this section. The crucial question is: what is the 
value of TLL parameters K v that enters into measured 
correlation functions? Obviously the answer cannot de- 
pend on the point where we arbitrarily decide to stop the 
RG flow, one has to keep in mind that gi terms usually 
are still finite at such a point. As discussed in Ref. [391] . 
the correlation functions arc of the form: 



r r/W 



exp 



y\\[i]di 



(29) 



dg 



2k r 



dl 



(2 - K pA - K aA )g 



2k F 



(28) 



It is the most relevant instability (already for K p a < 1) 
if one assumes SU(2) symmetric case in the spin sector 
when K a A = 1 , so naively we would expect it to dominate 
the low energy physics. However there are details that 
matter. First, the initial amplitude of this term is quite 
low. Second, the RG flow of this term will be strongly 
perturbed: 

• at high energies (A s» 0.2eU) the spin sector is still 
in the incoherent regime, certainly the dispersion 
is not yet linear (mind the value of J e ff)- Thus 
all terms which contain spin dynamics, like cos(0 CT ) 
will strongly fluctuate 

• the 2kp terms has to compete with Akp instabilities 
described in the previous section: large gz (in the 
initial part of the flow) and relevant g\ and g\ in 
the lower energies 



where \r\ = y 1 x 2 + {v u t) 2 is a distance in a time-space 
domain, W is an ultraviolet cut-off of the problem (~ 
2t), y\\[l] = K p — K p is a deviation of TLL parameter 
from the critical value K p of the considered flow. Thus 
the observed value could be interpreted not as K[l], but 
rather as weighted average of K[l] over longer and longer 
length-scales. In particular if we flow to a gapless phase 
then the observed K v corresponds to K*, which is then 
equal to the invariant of the flow A. 

Our problem is particularly difficult, because the an- 
swer to the above question strongly depends on yet un- 
known lowest energy physics. Indeed in our reasoning we 
have neglected several effects (t±, Dls{q = 2kp), prox- 
imity to Mott insulator, disorder) whose amplitudes are 
certainly smaller than A2 ~ 20meV and which seriously 
harm TLL below this energy scale. Based on this we 
can reassure the validity of our approach only down to 
- 20meU. 

Basically there are two distinct behavior which can oc- 
cur at these lowest energy scales. First, one of the above 
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mentioned perturbations can become pertinent just be- 
low the A2 energy scale and thus stop the RG flow. In 
fact all these perturbation in one way or another would 
not support further lowering of K p value, thus the fur- 
ther K p [l] can be then taken as vertical. The observed 
TLL parameter should then be close to if p [A2]. The sec- 
ond possibility, which is not unlikely, is that ID RG will 
remain valid down to much lower energies ~ lmeV (see 
Sec I VII Bj) . Then as indicated in FigJ3] the RG trajec- 
tory always stays close to separatrix. This statement is 
in fact quantified by the values of invariants of the intra- 
and inter-chain RG flows: A3 w A±, so it is always either 
<73 [/] or g±[l] keeping us close to separatrix. Then Eq|2T)l 
gives a power low with a K^. ho — 1/4 exponent times a 
logarithmic corrections in the form \og(r/W). 

We are able to compare our results only with the ex- 
perimental data corresponding to the 'high energy' range 
A G (0.02; 0.2)eV. In fact this is also the energy scale 
down to which, with no experimental doubt, the TLL 
persist. One should also keep in mind that the TLL a 
exponent should be computed within the ladder model 
introduced in Scc lVBI In the following we assume that 
it is only K p responsible for a 7^ 0, while K a = 1. Sev- 
eral different techniques have been used to measure the a 
exponent. These include PES and ARPES (4-150meV; 
5-300K), STM(0-50meV; 10-50K), resistivity (30-300K); 
where in parenthesis we have indicated temperature and 
frequency ranges where the fits were performed. We see 
that at least one of them is always larger than A2 which 
makes our theory applicable. All probes are consistent 
and gave the spectral function exponent a ~ 0.55 -7- 0.6 
which means K p 0.24 -j- 0.25. This is reasonably 
close to K p s,a[A-2] € (0.25,0.29) predicted theoretically 
in Sec lVTCl but also to K° p = 1/4. This shows the im- 
portance of the proximity to the Mott transition: it can 
cause a strong renormalization of K p TLL parameters 
towards lower values and thus observed large value of a 
exponent. 

One can be even more explicit: the only way to obtain 
values a > 1/2 (which translates into K p < 0.27) in our 
model is by assuming that the combined 4fcp perturba- 
tions (intra-chain g% together with g\ and g\ terms) can 
drive the RG flow and lead to strong renormalization of 
the charge TLL parameters. There is one more argument 
which supports this scenario. It is lack of spinon peak 
which has not been observed in any PES (or ARPES) 
experiment in the last two decades. This can be justi- 
fied theoretically only if a > 1/2, which provides strong 
limitations for the expected K* values. 

There is one experimental result which requires a com- 
ment. It is the temperature dependence of a invoked in 
Ref.[13], based on ARPES results. There the following 
value a = 0.9 was measured at high temperatures. This 
means (K pS + K pA )/2+{Kj + K-] L )/2 « 6 which trans- 
lates into K* ps 0.17. This value is much smaller than 
our predictions. We would like to emphasize that within 
TLL theory such behavior is very unlikely. To be pre- 
cise: as explained above it is not allowed to say that 



these different values comes out different points along 
RG trajectory. One model can possibly have only one 
ground state described by K*. In our opinion the ob- 
served temperature dependence can originate from a sig- 
nificant influence which phonons occupancy have on the 
LDA band broadening [17| The fact that a(T) dependence 
seems to change with sample preparation (see Fig. 3d in 
[20l |) strongly supports this interpretation. Then the real 
value of a is the low energy one. 

The other result is the estimate of a gap which accord- 
ing to our reasoning may be potentially opened by the 
AUf instabilities. However the value has been found to 
be extremely small ~ lmeV, thus our procedure, which 
relies on the RG analysis starting from 200meV, is insuf- 
ficient to make any definite claims. What is more this 
value is much smaller than A2, while there are several 
effects (like disorder or the interchain hopping tj_) which 
are potentially of order A2 and destructive for such a 
gap. If the hierarchy of amplitudes was opposite (gap 
much larger than t±) then the gap would suppress sin- 
gle particle hopping. In our case the outcome is unclear. 
Certainly a more sophisticated theory is necessary to un- 
derstand the unusual physics taking place below 20meV. 
Experimentally this is also a controversial issue. Such 
a gap (in the charge sector or in Atll{u)) has never 
been seen in any experiment probing uj > 3meV. So 
far the only exception are low temperature resistivity 
measurements Q, but even there a clear Arhenius-type 
activation behavior (at around lmeV) is preceded by an 
unusual power law like behavior. 



B. The validity of ID approximation 

It is known [2^] that the strong correlation effects (the 
formation of Luttinger liquid) are able to strongly reduce 
the value of t ± . To be more precise they strongly reduce 
the energy scale where system gains coherence along c- 
axis. With the value of the a exponent discussed above 
the perpendicular hopping tj_ , will get strongly renormal- 
ized down to an effective value: 

'i^O/*' ( 3 °) 

where £ = a + 1 is a single particle Greens function ex- 
ponent, and r/ t ± = 2 — £ is a scaling dimension of the 
hopping operator. This gives a suppression of t± by a 
factor ps 20. To be precise this means that due to the 
presence of strong interactions the hopping in the per- 
pendicular direction becomes coherent only ps lmeV. 

On the top of it there is also another source of t ± renor- 
malization, which originates from the competition with 
the terms V out (q = 2k F ) (as described in Sec |yiTJ|) . The 
2k p terms which are a functional of a charge asymmet- 
ric mode <f>pA (e.g. the one inducing a 7r— DW) tend to 
suppress t± [38(. Their influence is non-zero only when 
spin sector becomes coherent and t± is sufficiently small 
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(there is a Bessel function Jo(t±[l]) involved, which arises 
in a very similar way like the one in Eq fT7)) . It was shown 
[37l [H| that the significant suppression of tj_ may hap- 
pen only when t±[l] f=s V ou t(q = 2kp), which (according 
to scaling given in Eq |30l) can be the case in our problem 
but only for energies below A2. 

The frustration of the nearest and the next-nearest 
perpendicular hopping also suppress the coherence of 
perpendicular hopping also on second order when e.g. 
particle-hole processes are considered. This is immedi- 
ately visible if one writes a formula for any susceptibility 
x(q, u>) within the mean field, RPA level: 



for forward disorder is known in real space: 



A(x, t) = A T ll(x, t) exp 



D f K l 



(32) 



X FLL {q,u) 



1 + (2t± cos(^c/2) + 2f x cos(q ± c)) X TLL (q, u) 

(31) 



where we are using simplified version of perpendicular 
dispersion e(q~±) known from Eq lAll We see that the sign 
difference between t± and t' ± can cause a suppression 
of a second term in the denominator. This obviously 
weakens the q± dependence of susceptibility (thus also 
a corresponding observable), but also has its influence if 
one wished to develop a perturbative series to study the 
influence of t±. 

Finally there are also disorder effects which can lo- 
calize carriers in the perpendicular direction. They are 
described in the next section. 

From experiments (optical spectroscopy, Kadowaki- 
Woods ratio) we know that ID physics seems to be cor- 
rect even down to 2meV. In the light of the above dis- 
cussion these statements are not unreasonable. Thus, 
the assumption about validity of ID physics should be 
valid even down to few meV's and certainly in the " high 
energy" range of energies (200-20meV) this paper is dedi- 
cated for. The values of TLL parameters indeed strongly 
support the idea that the ID regime should be able to 
persist to temperatures much lower than the bare per- 
pendicular hopping. 



VIII. HOW PERTINENT IS SUBSTITUTIONAL 
DISORDER? 

The most likely source of a disorder in purple bronze 
are random vacancies of Li atoms. From the DFT 
results jl7j we know that energy shifts caused in this way 
within dispersive bands can be at most 15meV (when the 
Li are completely removed) . This sets the strength of the 
substitutional disorder potential. Because Li atoms are 
placed well outside zig-zag chains it is reasonable to as- 
sume the Coulomb potential interaction (~ q 2 ) between 
impurity and TLL, with predominant small q-exchange 
scattering. Thus the disorder will have primarily forward 
character and Df ss 15meV. 

Thus we can assume the model of TLL with a forward 
disorder to check if it explains the observed temperature 
anomalies of ARPES below 200K. The spectral function 



where Atll is the pure Luttinger liquid spectral func- 
tion, which is known e.g. see Eq|15l The behavior of the 
Fourier transform is easy to extract in the limits of large 
and small w. For the high-energy (15-150meV) range 
we expect the Atll power law scaling (because we work 
above the energy range disorder can affect). This is def- 
initely not seen in experimental data [!(| ■ For the low 
energies (2-15meV) we expect a convolution of standard 
TLL signal with an exponential decay. A broadening of 
Atll was indeed observed, but it was suggested[4lj that 
■the Gaussian convoluted with the LL Atll^ pro- 
vides better fit of the ARPES data. If one looks back 
how the eoi32lwas derived she/he finds that the approx- 
imation of uncorrelated scattering events is not obeyed. 
The presence of a Gaussian is like if the scattering events 
were not random, but momentum conserving. There is 
also discrepancy on the level of energy scales, experimen- 
tally the broadening has been observed already at 20meV 
which is larger than the maximal amplitude of disorder 
15meV. 

Let us now move to discussion of the lowest temper- 
atures effects. By analogy with a reasoning in Sec IV Al 
we expect that Df is accompanied with a component 
which provides source of random scattering events with 
large momentum exchange. Its amplitude should be one 
order of magnitude smaller and thus could potentially af- 
fect only the lowest energy scales. Taking into account 
the a exponent value (and repulsive character of inter- 
actions K p <C 1) the backward disorder, if present, has 
to be highly relevant. It is known that in ID systems 
the backward disorder, in the form of cosine operator, 
when present kills the superconductivity (of the BCS 
type) . Experimentally the unusual ARPES was observed 
precisely for the superconducting samples. As expected 
for more disordered samples the superconductivity dis- 
appeared, while the broadening of ARPES signal around 
100K was unchanged. The disorder also does not fit well 
neither with magnetoresistivity nor STM experiments. It 
is also claimed (from optical spectroscopy) that the mass 
of remaining mobile carriers decreases below 30K. 

In conclusion this simplest notion of disorder is incom- 
patible with observed effects. However the discussion is 
not yet closed. The strength of the disorder potential D f 
is comparable with the zero point motion or the thermal 
expansion effects, which will certainly affect its influence. 
The crystal structure is quite complex and supports much 
more sophisticated mechanisms e.g. the relative rota- 
tions of different Mo octahedra. These rotations provide 
another source of disorder in the material. However from 
the overall analysis of the structure it is certain that all 
these microscopic effects should affect t± stronger than 
the inside-chain Luttinger liquid physics. What is more 
t± « Df which means that relative strength of disorder 
is quite large. 
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IX. CONCLUSION 

We studied the low energy physics of a quasi- ID ma- 
terial, lithium molybdenum purple bronze. Already be- 
fore it has been shown that this material is extremely 
anisotropic, there is a linear dispersion along b-axis ex- 
tending down to 0.4eV, while in the perpendicular di- 
rection dispersion is at least two orders of magnitude 
weaker. In this study the physics is expressed in terms 
of field theory. Our work cover the energy range where 
the physics of charge modes, holons describes well the 
dynamics of the compound. It begins at around 0.2eV 
where the experiment has clearly shown an emergence of 
ID spectral properties, in particular the fermionic bands 
seem to merge at this energy scale giving rise to a single 
entity. The regime of our interest extends as low as ID 
physics is strictly valid. To be on the safe side we set it 
at around 20meV which is larger than any of the possible 
disturbances. 

The starting point of our study is the recent LDA-DFT 
result where the peculiar band structure of the material 
has been re-confirmed. Based on this we construct the 
effective low energy theory: the band structure around 
Fermi energy is casted into a tight-binding model. In 
addition a minimal model has to contain the strong cor- 
relation terms, laying beyond LDA-DFT approximation. 
The aim of the next few section is to provide the quanti- 
tative description of these strong interactions. 

We begin with parameterizing the interactions taking 
place inside a single zig-zag chain. From previous exper- 
imental and numerical works we are able to extract an 
effective real space model with a physically reasonable 
values of strong correlation parameters. Due to sparse 
arrangement of the chains the interactions have a finite 
range character, however in the first approximation we 
use U-V model to obtain the values of Luttinger liquid 
parameter K p . This estimate is based on several numer- 
ical works dedicated for models which are very similar to 
ours. All values converge at K p s» 1/3, which is rather 
low value, in particular it may allow 4fcp instabilities to 
dominate the physics. It also suggests that the umklapp 
processes are at play. In addition we also investigate the 
effects which can arise if one goes beyond simplistic U-V 
model. A separate section is dedicated to a spin sector, 
where we estimate basic energy scales. 

Later parts of this work are devoted to inter-chain 
physics. As the chains are grouped in a well distin- 
guishable pairs, it is tempting to propose a description 
within a ladder-like model. We use dielectric approxima- 
tion in order to estimate strength of inter-chain interac- 
tions, considering processes with both small and large 
momenta exchange. With this knowledge we propose 
that the description of physics in the considered charge 
regime should be done within framework of Luttinger 
liquid consisting of four modes, the two charge modes 
corresponds to the total (symmetric) and the transverse 
(asymmetric) fluctuations. This description fully incor- 
porates the inter-chain processes with small momentum 



exchange together with the intra-chain physics. 

The processes with large momentum exchange can be 
taken as perturbation and treated within RG approach. 
However, we claim that in order to achieve a valid de- 
scription of the system one has to consider the intra-chain 
umklapp together with them. We derive a full system of 
RG equations which cover intra- as well as inter-chain in- 
stabilities and study possible trajectory of the flow. This 
allows us to develop a description of purple bronze down 
to the limits of validity of ID theory. In the discussion we 
show that these limits can be safely extended to energy 
scales two or even three times smaller than 20meV. This 
allowed us to make a more extensive, quantitative com- 
parison with experimental results, which basically con- 
firms our theoretical insight into this complicated com- 
pound. 

We have achieved the effective low energy description 
of Lio gMoftOn compound which is able to explain con- 
vincingly experimental findings down to 20mcV. In par- 
ticular we showed that for the specific combination of pa- 
rameters, which are present in this material, an unusual 
situation may occur. Due to their mutual competition, 
the critical phase (the Luttinger liquid) is able to sur- 
vive down to the very low energy scales. Despite the two 
leg ladder formation, no gap opens in holons spectra (at 
least not above 20meV). This is contrary to the usual 
case where a significant gap is present in a transverse 
(asymmetric) mode of a ladder-like low-dimensional sys- 
tem. Such an unusual physics is in agreement with the 
physics extracted from experimental investigations. 

It is likely that for the lowest energies the purple bronze 
falls into category of doped Mott insulators with ex- 
tremely small gap, but further developments of the the- 
ory arc necessary to make any definite claims about this 
highly interesting regime. 



Appendix A: Derivation of tight-binding parameters 

The central result of DFT calculation is a band struc- 
ture of a given material. As it contains a huge amount 
of information usually it is difficult to deal with when 
the effective, low energy theory is constructed. In such a 
case a standard procedure is to approximate a solid by a 
tight binding model with a selected sites located at their 
positions Ti and unknown hopping parameters between 
them. Usually only the nearest t and the next-nearest 
t' neighbor hoppings are taken into account. We are in- 
terested in the low energy physics, thus the aim is to fit 
bands crossing Fermi energy Ep (in the given direction 
d) with the properly chosen parameters td t' d . 

As described in Sec|TT]in the case of Lio.gMoeOn most 
of low energy spectral weight in located on Afo(l) and 
Mo(4) sites (see Figflji)) so we take only them into fur- 
ther considerations. As explained only two bands cross 
Ep thus our aim is to fit these two dispersions. 

The dispersion relation for tight-binding model defined 
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above is known: 

e(k) = -tcos(k b b/2) - t' cos(k b b) 

- t±1 ^ t±2 cos(fc c c/2) - t±1 - t±2 sin(fc c c/2) 

- ^i±^g cos(fc c c) - ^~^ 2 sin(fc c c) (Al) 

where the first line describes the dispersion along b-axis, 
the last two the dispersion along c-axis (with the dimer- 
ization included). We neglected dispersion along a-axis. 
For the b-axis dispersion we assumed the simplest tight- 
binding model, as the LDA result suggest rather straight- 
forward interpretation of the bands. We introduced next 
nearest neighbor hopping t' mostly because the zig-zag 
chain structure seems to allow for this refinement, how- 
ever from the general shape of bands in b direction it does 
not seem to be necessary. 

For the c-axis dispersion the situation is quite different. 
We introduced much more terms because the curve is 
quite unusual, with a well pronounced double minima 
and a node at k c = 0, a feature quite difficult to fit within 
standard model. Based on the analysis of the structure 
presented in FigOJi) we deduce that: 

i there are two times more intra-ladder than inter- 

ladder links (they are linked either directly through 
Mo(l) — Mo(4) bond or auxiliaries through 
Mo(2) - Mo(5) bond) 

ii the inter-ladder hopping goes always though Mo{2) 

(or Mo(5)) atoms, these are two paths which can 
interfere. These hoppings are possible only every 
second site. 

iii the next-nearest neighbor hopping is allowed only 

through Mo(2) — Mo(5) thus should be much 
smaller than the one above. These hoppings are 
also possible only every second site. 

all of the hoppings goes either via a <5-bonds or 7r-bonds 
(the second are allowed only due to octahedra tilting). 
None of the links should be very much stronger than 
the others. From the crystal structure analysis all other 
hopping should be negligible, thus the task is to fit the 
dispersion with the above given model Eq lAll 

On Fig[6] we show an example of k c dispersion with 
conditions [i] — [iii] fulfilled. We see that it is possible to 
obtain quite similar shape to the LDA one provided that 
two interfering paths have the hopping parameters of the 
opposite sign. The only quantitative issue would be rel- 
atively large value of band splitting found in Ref . [l7| . 
It may appear from Hartree interaction term if the two 
bands had different orbital character. One should also 
remember that the value of splitting found in [Tsj was a 
bit smaller. We will leave this issue for further specialized 
studies like e.g. NMTO. 

These findings are summarized in a table Tab lllll 




FIG. 6: The dispersion along c-axis shown for varying direct 
S hopping along Mo(l) — Mo(4) bond. The strength of the 
other hopping path (through Mo(2) — Mo(5)) is taken as 0.02 
while the inter-ladder one 0.01. 



TABLE III: The values of parameters (given in meV) of the 
tight binding hamiltonian given in Eq lAll which fits best the 
LDA result [17J, the dispersion along c-axis. 



t±l 


t±2 


4: 




15 


5 


-10 


< 5 



Appendix B: Derivations of cosine terms 

We begin with a fermionic for of considered scattering 
processes: 1. the intra-chain umklapp: 

KosS = 53^(^^+92*1^-93^94*1 + 
9i 

^1^^4^+934^+94*1) ( B1 ) 

2. Interchain umklapp scattering (Akp ): 

HLs = 03" ^^+9llt^+922i^-932i^-94lt + 
9i 

^1^1^1^2^+932+^+94 It) (B2) 

where the following relation between momenta holds 
qi + qi — qi — qi = in order to preserve momentum con- 
servation during scattering; the low energy limit implies 
that for each momenta there exist the ultraviolet cut-off 

hi < A. 

3. Interchain exchange scattering (4fcp): 

H: os =gt kF f 'dxpf'pf* (B3) 
Jo 

4. Interchain exchange scattering (2k p): 

HZ os =gl kF f dx(%"f%» (B4) 
Jo 

where p 2 ^ F is the 2fc^ charge density in the i-th chain. 
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Usually one is also interested in correlation functions, 
which have the general form: 



Ro =< 2VO(x,t)O + (0,0) > 



(B5) 



where the simplest intra-chain examples of possible op- 
erators O are: 

• in Peierls channel - the charge density wave (CDW) 



C DW(X, t) = J2r,a,a> Vv(^, t) f <W 1pr(x, t) 

• in Cooper channel - the singles superconductivity 
(SS) 

SS {x,t) = Y, r ,a,a' ^r(x,t)S^l/) f (x,t) 

When we consider a single chain then a chiral fermion 
creation operator is related to bosonic spin and charge 
fields as follows (in the continuum limit): 



tpr.a(x) = }- — r? r ,g exp(irk F x) exp —1/ \/2(rcf> p (x) - 6 p {x) + a(r(f) a (x) - 6 a (x))) 
\2na 



(B6) 



where the coefficient r\ r ^ is the Majorana fermion which doesn't have any spatial dependance and it is introduced 
only in order to preserve anti-commutation for the fermion operators if). Usually they don't play any role in the 
physical description of the system but they are able to change signs of some correlation functions, so one has to take 
care about them. 

In the case of the ladder one has a straightforward generalization: 



i 

ipr,(r,v(x, t) ~ Vctm exp(irk F x) exp[- ~(r<f> +p + 9 +p + a(r<f) +(T + 6 +a ) + v(r4>_ p + 0_ p + <j(rcf)_ a + 9 -„)))] 



(B7) 



The above given equations Eq |B6l and Eq |B7l allow to 
rewrite all fermionic terms in the hamiltonian like e.g. 
the one in Eq lBll or in general any interesting operator 
into the language of bosonic fields. In our particular case 
we are interested in the following interaction terms: 

1. Umklapp scattering (at quarter filling) in the i-th 
chain 

Hcosa^gaf dxcos(2V8(f) ip (x)+5) (B8) 



where S indicates the doping away from a commensurate 
case, quarter filling in our case (where exp i(tt — Akp)x is 
not oscillating). 

2. Interchain umklapp scattering (4fcp): 

HLs = 9i [ dx cos(4cf>s P (x) + S) (B9) 



3. Interchain exchange scattering (4kp): 
HZ„ = qt kF I dxcos(A<f> Ap (x)) 



(BIO) 



4. Interchain exchange scattering (2k p): 



1 1 ~ 



Appendix C: Initial value of umklapp terms 

Estimating the value of <?3 at certain energy scale Ao 
is an impossible task, so all our results cannot be taken 
strictly quantitatively. This is a generic problem of the 
RG method present for any model even a simple half 
filled chain. The complexity of our system, makes the 
task even more tedious. 

We propose the approach, based on the first order ex- 
pansion of RG equations, to get a reasonable value of 
33 (Ao). However one has to keep in mind that as we are 
are working with rather large couplings, of order 0.1 (but 
always smaller than 0.25) so higher order terms can intro- 
duce non-negligible corrections even within single insta- 
bility flow. Thus the results has to be taken with caution, 
can be thought as only approximation. 

With this remarks being said we can proceed. One can 
in principle integrate out BKT equations to get a values 
of along the flow. One gets the following result [22l|: 



dx cos(2(f>A p (x)) cos(v / 2</>i -(a;)) cos(v / 2</>2cr(ir)) 



(Bll) 

These are the terms for which RG equations are de- 
rived in Sec lVIl In all the above we took a following con- 
vention: 

<f>i(x) ± (f>2 (x) 



9[l] 



y\\[i] 



A 



sinh 



Al + tanh _1 (A/j/jj) 



tanh 



Al + t&nh^iA/yV) 



(CI) 



(C2) 



V2 



(B12) 



where A is an invariant of the flow, in our case A = 1/12. 
Every point on an RG trajectory leads us to fixed point 
value K* ho . From numerics based on U-V model we 



which allows to go from one basis to another. 



found (Sec ED K* rho = 
quarter-filled chain K£ ho 



1/3 and we know that for 
= 1/4. Thus in our case 
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A = K* ho - K c rho = 1/12. Our aim is to extract the 
value of K r ho[l] — K^ ho + y\\[l] at certain energy scale 
Ai (corresponding to ps 0.2eV) which upon RG leads to 
known K* ho . This procedure can be thought as moving 
against the direction of the RG flow. 

The only missing quantity in Eq |C2l is yu[l = 0] which 
should be interpreted as the distance of an initial point 
of the flow K r ho[l = 0] (formally in the infinite energy, 
i.e. somewhere around UV-cut-off of the model) from 
K° ho = 1/4. Although around the fixed point the flow is 
rather vertical, we decided to take the most conservative 
(and giving the most modest value of g% (Ao)) assumption 
that K rho [l = 0] = 1/2 => y\i [I = 0] = 1/4. It is because 
certainly at K r ho = 1/2 the quarter filled chain remains 
insensitive to Mott localization. 

With this we are able to estimate the bare 173 (Ai) A to 
be around 0.35ey. Then yii(Ai) is only 6% larger than 
corresponding g(Ai), while the position of the separa- 
trix of RG flow in the first order is set by y\\ — g. The 
distance from the separatrix is of the same order like ad- 
ditional intra-chain terms described in Sec lIV Cl From 
the remarks above we know this is only an estimate, but 
it certainly implies that the umklapp, even in intra-chain 
U-V approximation, is an order of magnitude larger than 
any of V±(q = ikjr). Although it seems to be irrelevant, 



but its bare amplitude is large and quite close to sepa- 
ratrix so the intra-chain umklapp has to be taken with 
care in the RG analysis below. 

The fact that bare 173 falls very close to separatrix can 
be confirmed by an independent reasoning. The posi- 
tion of separatrix is frequently constrained by symmetry 
considerations. Taking into account the hierarchy of en- 
ergy scales, we take simplistic approximation of spinless 
fermions in U — > 00 (like in Sec lTVB|) . This obviously 
overestimates 53, but numerics [29( shows us that be- 
yond a threshold U = At, the dependence on U is weak, 
thus the corrections breaking (particle-hole) symmetry 
should be quite small. The advantage is that now one 
can map charge sector on pseudo-spin model (empty and 
occupied sites), for which the position of separatrix of 
RG flow (phase transition) in terms of bare parameters 
is known exactly thanks to SU(2) symmetry. It is located 

at V?J[* = 2t e ^ where V^[* accounts for all unscreened 
interactions along the chain (while t e ^ w t, if one want 
to be more precise he would found it slightly reduced with 
respect to t because of dimerization). If we put values 
found before (in table Tab HH for V*f f ) we indeed find 
ourself very close to ■ 
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